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B. J. MacLennan 
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Monterey, CA 93943 

Abstract: 



In this report we demonstrate the relational programming language RPL by using it to develop 
four programs. These programs are: (1) computing a table of word frequencies from a text; (2) 
minimizing a deterministic finite state automata; (3) Gaussian elimination; and (4) a simple data 
processing example involving updating an employee file. Appendix A shows transcripts of execu- 
tions of the programs on the Browm and Mitton interpreter Brown&Mitton . The reader is 
presumed to be familiar with RPL, which is described in MacLennan83 . For convenience, how- 
ever. Appendix B contains the RPL grammar, and Appendix C describes the language accepted 
by the Browm and Mitton interpreter. 



1. Computing Word Frequencies 



The first example, which is adapted from !MacLennan83 , it to compute a frequence table F from 
a text (sequence of words) 5. That is, given 5 such that S[i is the i th word, we compute F such 
that F[w is the frequency (number of occurrences) of word w in 5. For an example, we take 5 = 



<“to”. “be’\ “or”, “not”, “to’ 



5 = 



be” >, w r hich is just an abbreviation for the relation: 
{ l:“to”, 2:“be”, 3:“or”, 4:“not”, 5:“to”. 6:“be”} 



Pictorially, 



1 


“to” 


2 


“be” 


3 


“or” 


4 


“not” 


5 


“to” 


6 


“be” 



In this case, the desired frequency table F is: 

F = {“to”:2, “be”:2, “or”:l, “not”:l} 

Pictorially, 



“to” 


2 


“be” 


2 


“or” 


1 


“not” 


1 



In other words, “to” occurs twice, “be” occurs twice, “or” occurs once, and “not” occurs once. 
Of course, since F is a relation, the order in which the elements are listed is irrelevant. 

To develop the general word-counting program, we work through this particular example. 
Since the frequency table is a function from the words to their frequencies, the first step is to 
reverse columns of S: 



7-1 - 



{“to”:l. “be”:2, “or”:3, “not”:4. “to”:5. “be”:6} 



* The work reported herein was supported by Contract N00014-85- WR-24057 from the Office of Naval Research. 
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Pictorially, 



“to’" 


1 


“be’' 


2 


“or” 


3 


“not” 


4 


“to” 


5 


“be” 


6 



Notice that this relation is not a function (i.e., it is not single valued). We can make it a func- 
tion by forming the ‘‘unit image” of the table: 1 



unimage S 1 = {“to”:{l, 5}, u be”:{2, 6}, u or”:{3}, u not”:{4}} 



Pictorially, 



“to” 


{1.5} 


“be” 


{2.6} 


“or” 


{3} 


“not” 


{4} 



This tells us, for example, that the word u be" occurs in positions 2 and 6 in the text. We do not 
need to know the places where a given word occurs, but only the number of such places. There- 
fore, we send the preceding table through the size (cardinality) function (by the relative product' 
operation): 

unimage size = { u to”:2, “be”:2. “or”:l, u not”:l} 

This is the desired result; the final step is depicted in Figure 1. Notice that since size is defined 
for all sets, it is in effect an infinite relation: this is permitted in RPL. 

The resulting program is: 

F = unimage S~ l size 

We can turn it into a function definition to compute the frequency table for any text 5 by: 

freq S = unimage S~ l size 

It remains to define the ‘unimage 1 function, which is not built into RPL. On the other hand, 
RPL does have the builtin operator unimg, defined to that T unimg x is the set of all y such 
that x:y E T. This can be used to define unimage. To see this, note that the left section 
T unimg is the function that takes any x into its image under T. Although T is finite (and 
extensional), sections are always intensional, so it is necessary to to convert T unimg to its 
extensional equivalent. This is accomplished with the RPL restrict operation, which converts an 
intensional relation to an extensional relation by restricting its domain to a finite set. Hence we 
define: 

unimage T = dom T restrict T unimg 

The following is an example RPL session that defines the freq function an applies it to a particu- 
lar text (‘?>’ is the RPL prompt): 



1. This function is related to the RPL unimg operator: see below. 
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“to” 


{1. 5} 


“be” 


{2.6} 


“or” 


{3} 


“not” 


{4} 



0 


0 


{1} 


1 


{4} 


1 


{2,6} 


2 


{ 4 , 8 } 


2 


{ 4 , 6 , 7 } 


3 



“to” 


2 


“be” 


2 


“or” 


1 


“not” 


1 



Figure 1 . Piping unimage S 1 Through size Function 



?> unimage T = dom T restrict T unimg 
?> freq S = unimage S~ l size 
?> freq <“to”, “be”, “or”, “not”, “to”, “be”> 
{“be”:2, “to”:2, “not”:l, “or”:l} 



?> done 

The appendix contains the actual transcript of this RPL session; it shows how programs must be 
represented for the Brown and Mitton interpreter. The preceding formulas were produced from 
this transcript by a pretty printer. 

2. Minimization of Deterministic Finite Automata 

The next example program is the equivalence and minimization of deterministic finite automata 
by an algorithm developed by Robert Floyd 2 . We assume that we have a finite alphabet E and a 
finite set Q of states. The set F C Q represents the final (accepting) states. The finite relation 
T is such that for aGE, T^a is the transition relation for the symbol a. That is, < <7, q' > E T 4 a 
if and only if the symbol a takes state q into state q ' . 



2. Private communication, 1985. 
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Our goal is to define a relation R ^ such that << 7 , q > £ R x if and only if q and q' are not 
equivalent states. This is done in a series of steps, starting from pairs of states that are known to 
be inequivalent, namely the final and nonfinal states: 

R 0 = Fx(Q\F) 

We now work backward: any states that under the same input lead to inequivalent states are 
themselves considered inequivalent. For example, R 0 relates inequivalent states; R x relates states 
related by R 0 together with those that under the same input character are taken into states 
related by R 0 ; Rn relates states related by R 0 together with those that under the same one or two 
input characters lead to states related by R 0 \ and so on. Each step of this process is accom- 
plished by a function xp: that is, we will define xp so that 7 ? |+1 — xp R x . It will be easy to see that 
this process converges in n = (size Q)~ steps, so 

Roc = r R 0 

Next we consider xp. 

As a preliminary we define the polymorphic image of one relation under another relation. If R 
and 5 are two relations, then the polymorphic image under R of 5. or more briefly the R- image 
of S, is defined 

R S = R S j R~ l 

This has the following property: R ± S relates x to y if and only if there are u and v such that R 
relates 1 to u and y to t\ and 5 relates u to v. That is, Cr, y> £ R ± 5 if and only if there are u 
and v such that <x. u> 6 R. <y. v> E R and <t/, v> £ 5. This can be visualized: 



R 

1 ^ u 



R 



S : 

V 



y 



TT* 



s 

v 



Now, if we have that R t relates inequivalent states, then (T * a) _L R t will relate those states 
that are carried by symbol a into states inequivalent by R { . Thus, we define /? l+1 so that it 
relates those states that are related by R x together with those related by (T [ a) ± /?,. for any 
Now% if by [j_ /?, w'e mean the function that takes a polymorphic image of R^ that is, 

l-L *,]* = * ± R> 



then it is easy to see that 

m LR,])la = [iRMTia) = (T^a)±R, 

Hence, the union of ( T j a) _ R r for all a E E, is just the union of the range of the relation 
T t j_ R t . This yields the definition of R t and hence xp: 

R t+ i = ^ R, = R< u u( rn g ( t ! i-L ^,])) 

This completes the definition of the inequivalence relation R 00 . Two states are now' equivalent if 
they are not inequivalent: 

R = = Q 2 \(R x DRp) 

where Q 2 means Q X Q. 



The minimal machine is constructed on the basis of the equivalence classes of states under 
R^. The equivalence class of a state q is just the unit image under of q: 



- 4 - 



eclass q = R = unimg q 



The set of all such equivalence classes results from taking the image of Q under eclass: 

Q = = eclass img Q 



In general, we define 



equiv = eclass img . 

so Q = = equiv Q. We take Q = to be the states in the minimal machine. 

It remains to construct the transition relation T = of the minimal machine. 
<q. q > E T l a we want 



For all 



< eclass q , eclass q' > G T _ ^ a 
Thus T = l a is the isomorphic image under eclass of T j a: 

T = l a = eclass $ ( T [ a) 

Thus is the (finite) composition of eclass $ and 7\ which is the (finite) relative product of T 
and eclass S : 

T = = T [eclass $] 

The remainder of the minimal machine is easy to construct. For example, the final states are just 
the equivalence classes of the original final states: 

F = = equiv F 

There follows the actual relational program to minimize a small automaton. It makes use of two 
auxiliary functions a and p for defining the union of a set of sets: 



— DFA Minimization 



— Utility Functions 

1st = \l 1] 

2nd = [l 2] 
r l s = r | s | r -1 

O f = {f 0 (lst7(f 0 2nd)))7([\] 0 (17 (un » f)) » 2nd)) 

/ p i = 1st • (c f while ( 0] • 2nd)) ° i , 

u = M p 0 

-- Example DFA 

S = {1,2} 

T = {1 : {10 : 10, 20 : 20}, 2 : {10 : 30, 20 : 30}} 

Q = {10,20,30} 

F = {30} 



— Minimization 



Q- 


= 


Q x Q 


n 


= 


size Q- 


R 0 


= 


F x (Q \ F) 


v R 


= 


R u( rn g ( T i i-L R\)) 


Roc 


= 


ll> n Rq 


R = 


= 


Q 2 \(R x uR-J) 


eclass 


= 


R = unimg 


equiv 




eclass img 


Q- 


= 


equiv Q 


T = 


= 


T | eclass S 


F m 


= 


equiv F 



— Minimized DFA 
val <? = 

{ {10. 20}. {30} } 
val T - 

{ 1 : { { 10 . 20 } : { 10 . 20 } }. 

2 : { {10. 20} : {30} } } 
val F = 

{ { 30 } } 

The val command prints the value of an identifier. 

3. Gaussian Elimination 

The matrix is represented as a vector of vectors: 



<< 


a n - • • 


• ’ fl ln- 




>, 


< 


floj * . . 




b. 


>, 


< 


a nl > • ■ 


’ * T 




>> 



For the sake of simplicity we assume all the a j; are nonzero. We use the operator 4’ to select the 
kth element of a vector. Thus "M j k ’ is the k th row of M and | [j k] J is the k th column of 
A/. 

The Gauss Elimination function will use n successive steps. Each of these steps will accom- 
plish the transformation 

<A/. k> -+ <A/\ k + \> 

where M ' is obtained from M by performing the elimination process on the k th column: 

M' = elim <A/. k> 

Thus the complete process is defined: 

Gauss M = (elim for <1 , ... , n>) M 

Here we make use of the functional l (/ for S) x' which computes the sequence of values 

y, = /< i. S x > 
y 2 = /< j/p S 2 > 

Vn /' — 1 ’ 
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and returns y n . The Tor' functional is defined in terms of reduction as follows: 

/ for 5 = l@ 5 • [/ § 

This can be understood by the expansion: 

(/ for 5) x = ([@ 5; • !/ §]) x = (@ Sj (,/ §] x) 

= (@ 5] (/ § x) = (/ § x) @ 5 
= (/§*)$ 

We turn now to the elimination process. We want Al' = elim <M, k> , where M' results from 
M by zeroing all entries in column k of A/, except the entry in row k , which is set to one. This 
can be accomplished by subtracting an appropriate matrix E from M: 

M' = mat d if <A/. E> 



Here ‘matdif is a component-wise matrix difference function. 



The matrix E is produced by multiplying the appropriate factors by the individual rows of M. 
For the first elimination step E is: 

< (1 — 1/ a^) A/j, 

( a 2\ ! a ll)^l’ 

E i = : 

( fl n i/ a n)^i > 

The matrix resulting from subtracting E x from M is 



<< 


1. 


fl JO , 


'• a ln'- 


b l >, 


< 


0 , 


/ 

floo , 


a 2n 1 


b 2 1 


< 


0, 


a n2 ' 


a nn ' 


b \ » 



At the next stage the elimination matrix is: 

< ( a 12 / a 22 2 i 

(1-1/ a 22 ') A/o \ 

e 2 = : 

( a n 2 / a 22 )^2 > 

In general, if M' ' is the matrix resulting from the k — 1st elimination step, then the elimination 
matrix for the k th step is 



E k “ (l a *A: l}/ a kk )Mk ’ 



i a nk ' ' / a kk ' ')M k > 

It is easy to see that E k results from multiplying a vector V k by the A: th row of M' ' . This is just 
the outer product of and the k th row of M' ' : 

E k = outerprod <V k * M' \k> 



The vector V k is 



< 



L \k 



/ a kk 



^ k ~ i a kk l)/ a kk • 



a nk"/ a kk" > 

This is obtained by forming the scalar product of 1 / a kk ' ' and the vector 

k* = <a \k * a 2k y • • • * a kk “ 1 * • * ? a n k > 

This in turn is the result of subtracting from the A: th column of \1' ' the unit vector unit <A/, £>, 
which has a 1 in the k th position, and a 0 in all others. 

We now develop an explicit relational formula for V k . For this purpose it will be convenient 
to treat it as a binary function V k = V <A/, k>. We have: 

V <A/, k> = scaprod <1 / diag <M, > . U k > 

= scaprod <1 diag <A/, k>. vecdif <column < A/, &>, unit <A/, k>>> 

The parameter <A/, A:> can be factored out by use of the construction operation ‘T\ which is 
defined so that [f~g)x = <fx. gx>. Factoring, we have: 

T <A/. k> = scaprod ^1 / diag <A/, /:>. vecdif <column <A/. /:>, unit <A/. k>>> 

= scaprod <(jl /] 0 diag) <A/, /:>. vecdif ((columnTunit) <A/. /: >)> 

= scaprod <(|l /I 0 diag) <Af, /:>, (vecdif 0 (column , unit)) <A/, /:>> 

= scaprod ((( 1 / 0 diag) , (vecdif 0 (column , unit))) < A/, £>) 

= (scaprod ° ((1 / ° diag)T(vecdif ° (column“Tunit))) <A/. k> 

Canceling <A/, k> from both sides yields an explicit formula for V: 

V = scaprod 0 ((:1 /] ° diag)T(vecdif 0 (columnTunit))) 

We proceed similarly to get a formula for elim: 

elim < M, /:> = matdif <A/, £’ jt > 

= matdif <A/, outerprod < M [ k>> 

= matdif <A/, outerprod < V < A/, /:>, A/ | /:>> 

We perform some minor rearrangements so that <A/, k> can be factored out of the right-hand 
side: 

elim <A/. k> = matdif <A/, outerprod < V <A/, /:>, [|] <M , k>>> 

= matdif <M, (outerprod 0 ( l~[i!)) <Af. A>> 

= matdif <[j 1 <A/, /:>, (outerprod ° ( VT[1])) <A/, /:>> 

= (matdif « ((| ljT(outerprod 0 0"7[i])))) <M^ k> 

Canceling <M , A> from both sides yields an explicit formula for elim: 

elim = matdif ° ([{ l]T( ou terprod ° ( KT[ij))) 

A complete RPL session demonstrating the Gaussian elimination function follows: 
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— Utility Functions 

con k = 
transmap / = 
vecdif = 
scaprod <k, v> = 
outerprod <u, v> = 
matdif = 
column <A/, k> = 
unit <A /. k> = 
diag <A/. /:> = 

/ for S = 



X x k 

[I /] • \f. 

transmap — 
v k x ] 

u (scaprod 0 , v\) 

transmap vecdif 

M [i k] 

< 1 , . . . . size M> 
M l k l k 

[@ • [/ §■ 



[[=*] — con 1: con 0 



— Gaussian Elimination 



V = scaprod ° (( 1.0 /] c diag)T(vecdif ° (columnTunit))) 
elim = matdif 0 ( ^ 1 T(outerprod ° ( ^T[i] ))) 

Gauss M = (elim for < 1 , . . . , size A/>) M 

— Example Matrix 

M = << 3 . 9 , 33 >. < 2 , — 1 , 1 >> 

— Execution 



Gauss M 



<< 1.0, -2.38419E-7, 2.0> 

< 0.0, 1.0, 3.0>> 

The matrix A/ represents the equations 

3x 4 - 9y = 33 
2x - y = 1 

The result of Gauss M correctly reflects the solution x = 2, y = 3. 



4. Employee File Update 

Next we consider a simple data processing example adapted from MacLennan83j. We are given 
an employee file F indexed by employee number. That is, F [ n is the record for employee 
number n. The employee records themselves are represented by functions from attribute names 
into attribute values. For example, if R is an employee record, then R j “N” is the employee’s 
name, R j U R” is his hourly rate, and R j “H v is the hours worked this pay period. Here is an 
example employee file containing three records: 

F = {124 : {“N” : “John”, “R” : 10, “H” : 100}, 

118 : {“N” : “Bill”, “R” : 15. “H” : 120}, 

207 : {“N” : “Sally”, “R” : 14, “H” : 115}} 

We are also given an update file U such that U [ n is the number of hours worked this week bv 
employee number n. For example: 

U = {118 : 6, 124 : 40, 207 : 40} 

Our task is to generate an updated employee file F' in which the hours worked (“FI”) field has 
been updated. 

First we define ‘sumhrs’ so that if R is an employee record and h is the hours worked this 
week, then sumhrs </?, h> is the new total hours. Clearly, 
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sumhrs </?. h> = (R | "H”) + h 

Alternately, we can define this function variable-free style: 

sumhrs = [+] 0 ([„ "H”] I) 

It is easy to see the two are equivalent: 

sumhrs </?, h> = ([+] - ([l “H”] |j I)) </?, h> 

= [+] (([i “H”] || I) </?, h>) 

= [+] <[1 “H v ] R , I h> 

= (R l “H”) + h 

Our next task is to replace the old value of “H” field by h' . This can be accomplished by the 
ordered union operation V* For example. 

{“H” :h'};R 

will return a record R' in which R' j "H” = h' but all other fields of R' are the same as in R. 
How do we get the relation {"H*" : A'}? Since this is just a sequence that’s equivalent to the 
array <"H’\ h' >. we can use as to convert the array to a sequence. We solve for the function / 
that computes {'*H” : h'} from </?. h> as follow's: 

/ </?, h> = {"H” : sumhrs </?. h > } 

= as <“H V . sumhrs <R. h>> 

= as ( “H” J (sumhrs </?, h>)) 

= (as ° ' U H’’ ,j 0 sumhrs) </?, h> 

Hence, 

/ = as ° [ U H‘* ,J ° sumhrs 

It’s necessary to get the corresponding records from the F and U files together so that they can 
be processed by /. This is accomplished by the extensional construction operation # defined so 
that (F f U) ; n = <F ^ n. U i n>. With the given example files we have: 

F # U = { 124 : <{“N” : "John'’, “R” : 10. "H ,? : 100}, 40>. 

118 : <{"N V : "Bill’*, U R” : 15, "H” : 120}, 6>. 

207 : <{“N’’ : "Sally”. U R” : 14, U H” : 115}, 40>} 

Notice that the pairs <F j n. U | n> are just the inputs required for /. We combine the preced- 
ing results into a update file ‘upd’ defined so that upd j n is {“H” : /T}, representing the new' 

hours worked for employee number n. In this case, 

upd = { 124 : {“H” : 140}, 

118 : {"H” : 126}, 

207 : {"H” : 155}} 

It’s easy to solve for upd by using the relative product: 

upd ■ n = / ((F # U) i n) 

= {(F # U) ! /) : R 

Hence, upd = (F f U) \ f . Substituting for / yields: 

upd = (F f U) (as ° ["H” ,] ° sumhrs) 

Now' we’re almost done. W’e w^ant each record in F' to be the ordered union of the corresponding 
update record in upd and old record in F. Hence we solve: 
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F' ± n = (upd * n) ; (F j n) 

= [;] <upd „ n, F n> 

= !:] ((upd f F) ' n) 

= ((upd # F) | [;:) . n 

Hence, 

F' = (upd # F) [;, 

The complete session follows: 



— The Files 



F = { 


124 : 


{“N v : 


"John" 


. “R” 


: 10. 


“H” : 100}, 




118 : 


{“N” : 


“Bill”, 


“R” : 


15, “ 


H” : 120}, 




207 : 


{“N” : 


“Sally’ 


“R" 


: 14. 


“H” : 115}} 


U = { 


118 : 


6. 124 : 


40. 201 


' : 40} 






Computing the New File 








sumhrs 


s [ 


+ ] * (U 


“H”] |l 


I) 






upd 


= ( F # U) 


| (as ° 


“H” 


, ° sumhrs) 


F' 


= ( 


upd # F) | [;] 








The New 


File 












val F' 














{124 : 


{“FT 


’ : 140. 


“N” : “ 


John” 


. “R‘ 


’ : 10}, 


118 : 


{“H - 


' : 126. 


“N” : “ 


Bill”, 


UR” 


: 15}, 


207 : 


{‘•H - 


' : 155. 


“N” : “ 


Sally” 


, “R’ 


14} } 



This result correctly reflects the fact that John (employee 124) has worked 124 hours. Bill 
(employee 118) has worked 15 hours, and Sally (employee 207) has worked 14 hours. 

It is simple to modify the program so that it uses the input files OldMaster and Updates, and 
defines the output file NewMaster: 

F = file “OldMaster” 

U = file “Updates” 
sumhrs = (+] • ([„ U H”] | I) 

upd = (F f U) (as e “H” , ° sumhrs 
file “NewMaster” = (upd j* F) I [; 

5. References 

Brown&Mitton Brown, J. R., and Mitton, S. J., Relational Programming : Design and Imple- 
mentation of a Prototype Interpreter , MS thesis, Naval Postgraduate School, June 1985. 

MacLennan83 MacLennan, B. J., “Relational Programming,” Naval Postgraduate School Com- 
puter Science Department Technical Report NPS52-83-012, September 1983. 
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APPENDIX A: EXAMPLE RPL SESSIONS 

This appendix contains transcripts of actual RPL sessions with the Brown and Mitton inter- 
preter. Note that the interpreter follows the Interlisp convention of permitting a bracket *]’ to 
close any number of open parentheses. 

Example 1: Word Frequence 

It will be seen that the RPL interpreter computes a relation containing redundant tuples. They 
do no harm, but can be eliminated (by a quadratic algorithm) if desired. The transcript follows: 

Loading RPL DO YOU WANT TO RESUME A PREVIOUS RPL SESSION? <y/n> n 

RPL INTERPRETER ON LINE!! 

?> S == (list "to" "be" "or" "not" "to" "be"] 

?> (S sup -1] 

( rel ( be 6 ) ( to 5 ) ( not 4 ) 

( or 3 ) ( be 2 ) ( to 1 ) ) 

?> unimage T == ((dom T) restrict (Isec T unimg] 

?> (unimage (S sup -1] 

(set (be (set 6 2 ) ) (to (set 5 1 ) ) 

( not (set 4 ) ) (or (set 3 ) ) 

(be (set 6 2 ) ) (to (set 5 1 ) ) ) 

?> ((unimage (S sup —1)) rp size] 

( rel ( be 2 ) ( to 2 ) ( not 1 ) (or 1 ) 

( be 2 ) ( to 2 ) ) 

?> freq S == ((unimage (S sup —1)) rp size] 

?> (freq (list "to" "be" "or" "not" "to" "be"] 

(rel (be 2 ) (to 2 ) ( not 1 ) (or 1 ) 

( be 2 ) ( to 2 ) ) 

?> done 

DO YOU WANT TO SAVE ENVIRONMENT FOR FUTURE USE? <y/n> n 
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Example 2: Minimizing DFA 

For this example we assume that commands for defining the DFA and performing the minimiza- 
tion are on a file, u examples/dfa.rpl’\ whose contents are: 

(1st == (rsec sel 1)) 

(2nd == (rsec sel 2)) 

(r ppd s == (r | (s | (cnv r)))) 

(sigma f == ((f o (1st (. bar) (epsilon o 2nd))) (, bar) ((op ) o ((I (. bar) (un o epsilon)) o 2nd)))) 

(f rho i == (1st o (((sigma f) while ((rsec != empty) o 2nd)) o (Isec i ,)))) 

(union == ((op cup) rho empty)) 

(SIGMA == (set 1 2)) 

(T == (rel (1 : (rel (10 : 10) (20 : 20))) (2 : (rel (10 : 30) (20 : 30))))) 

(Q == (set 10 20 30)) 

(F == (set 30)) 

(Q sup_2 == (Q cart Q)) 

(n == (size Q sup_2)) 

(R_sub_0 == (F cart (Q F))) 

(psi R == (R cup (union (rng (T rp (rsec ppd R)))))) 

(R sub_ inf == ((psi sup n) R sub_0)) 

(R sub = == ((Q sup_2 R sub__ inf) cap (Q sup__2 (cnv R_sub inf)))) 

(rom eclass =— (Isec R s u b = unimg)) 

(rom_equiv == (Isec rom_eclass img)) 

(Q sub = == (rom equiv Q)) 

(T sub ,= == (T rp (Isec rom eclass $))) 

(F sub = == (rom equiv F) ) 

EOF 

This file is executed by being loaded into RPL. The resulting transition function and states of 
the minimal machine are then displayed. They can be seen to be sets of sets, since the states in 
the minimal machine are represented by equivalence classes. 0 The transcript follows: 

DO YOU WANT TO RESUME A PREVIOUS RPL SESSION? <y/n> y 

INPUT FILENAME 

examples/dfa.rpl 

Loading Session loaded 

?> val Q sub = 

( set ( set 10 20 ) ( set 10 20 ) ( set 30 ) ) 

?> val T sub ,= 

( rel ( 1 ( rel ( ( set 10 20 ) ( set 10 20 ) ) ( ( set 10 20 ) ( 
set 10 20 ) ) ) ) ( 2 (rel ( ( set 10 20 ) (set 30 ) ) ( ( set 10 
20 ) (set 30 ) ) ) ) ) 

?> val F sub = 

( set ( set 30 ) ) 

?> done 



3. Note that as usual there is benign redundancy in the sets. 
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Example 3: Gaussian Elimination 

The program for performing the Gaussian elimination is in the file “examples/gauss.rpl v , whose 
contents are: 

(con k == (func x k)) 

(transmap f == ((rsec rp f) o (op #))) 

(vecdif == (transmap (op — ))) 

(scaprod (k v) == (v rp (Isec k times))) 

(outerprod (u v) == (u rp (scaprod o (rsec . v)))) 

(matdif == (transmap vecdif)) 

(column (M k) == (M rp (rsec sel k))) 

(unit (M k) == ((listrange 1 to (size M)) rp (if (rsec = k) -> (con 1) ; (con 0)))) 

(diag (M k) == ((M sel k) sel k)) 

(f for S == ((rsec @ S) o (Isec f red))) 

(V == (scaprod o (((Isec 1.0 divide) o diag) (, bar) (vecdif o (column (, bar) unit))))) 

(elim == (matdif o ((rsec sel 1) (, bar) (outerprod o (V (, bar) (op sel) ))))) 

(Gauss M == ((elim for (listrange 1 to (size M))) M)) 

(M == (list (list 3 9 33) (list 2 -1 1))) 

(a == (diag (list M 1))) 

(b == (vecdif (list (column (list M 1)) (unit (list M 1))))) 

(v == (scaprod (list 0.33 b))) 

EOF 

The session shown in the following transcript performs the Gaussian elimination on the matrix 
M: 

DO YOU WANT TO RESUME A PREVIOUS RPL SESSION? <y/n> y 

INPUT FILENAME 

examples/gauss.rpl 

Loading Session loaded 

?> (Gauss M] 

( rel ( 1 ( rel ( 1 1.0 ) ( 2 -2.38419E-07 ) ( 3 2.0 ))) 

( 2 ( rel ( 1 0.0 ) ( 2 1.0 ) ( 3 3.0 ) ) ) ) 

?> done 

Note that the resulting matrix is printed as a relation rather than a list of lists, since it is quite 
expensive for the interpreter to determine if a relation is in fact a list. 
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Example 4: Data Processing 

In this example, the employee file to be updated is small (thrae records), and so typed in interac- 
tively. More typically, the RPL file facility would be used to load F from disk. The transcript 
follows: 

RPL INTERPRETER ON LINE!! 

?> F == (rel (124 : 

(rel ("N" : "John") ("R" : 10) ("H" : 100))) 

(118 : 

(rel ("N" : "Bill") ("R" : 15) ("H" : 120))) 

(207 : 

(rel ("N" : "Sally") ("R" : 14) ("H" : 115] 

?> U == (rel (118 : 6) (124 : 40) (207 : 40] 

?> (F # U] 

( rel ( 124 ( rel 

(1 (rel ( N John ) ( R 10 ) ( H 100 ) ) ) 

( 2 40 ) ) ) 

( 118 ( rel 

(1 (rel ( N Bill ) ( R 15 ) ( H 120 )) ) 

(26))) 

( 207 ( rel 

(1 (rel ( N Sally ) ( R 14 ) ( H 115 ) ) ) 

( 2 40 ) ) ) ) 

?> sumhrs == ((op +) o ((rsec sel "H") || I] 

?> upd == ((F # U) rp 

(as o ((Isec "H" ,) o sumhrs] 

?> F - == ((upd # F) rp (op ;] 

?> val F' 

( rel ( 124 

(rel ( H 140 ) ( N John ) ( R 10 ) ) ) 

( 118 

(rel ( H 126 ) ( N Bill ) ( R 15 ) ) ) 

( 207 

(rel ( H 155 ) ( N Sally ) ( R 14 ) ) ) ) 

?> done 
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APPENDIX B: RPL GRAMMAR 



session = command done 

prefixid [ identifier ] = expression I 



command = 



expression 



application = 



display expression 

I expression infix application \ 
superscription I 

( application primary 

iter ' ' primary -> primary 

application 



superscription = expression sup 



literal 

prefixid 



primary — 



infix 

infix primary 
primary infix 

primary primary ; primary 

( expression .. expression ) 

{ expression .. expression } 

< primary , • * • > 
file string 



infix = infixop bar 
identifier = letter 



letter 

digit 



prime 



prime = 



literal = < 



digit * . digit ^ E [ -f- | — ] digit 
string 

true 

false 



string = u char ’* 
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infixop = 

sel . : cup member nomem [subset subset = -> <- restr : cl cr cap \ 
@hat ! cat @ . | $ red + - times divide ! = < > <= > = 

andsign orsign cart 



prefixid 



identifier 

prefixop 



prefixop = 

- un cur unc theta size str DELTA inv dom rng mem Lm Rm Mm run lun bun 
ini t term alpha omega ALPHA OMEGA min max mu index select join as sa saO 
rp rpi rsort sort unimg all ssm img curry uncurry PHI Id while upsilon 
phi delta PI extend restrict wig not 
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APPENDIX C: RPL INPUT FORM SUMMARY 



TABLE 1 . Primitive Extensional Operations 



Name 


Old Input Form 


New Input Form 


Publication Form 


selection 


t sel x 


t sel x 


t i i 


relative product 


t u 


t | u 


t | u 


construction 


t . bar u 


t # u 


t # u 


pair formation 


x : y 


x : y 


x : y 


union 


t cup u 


t cup u 


t U u 


unit set 


un x 


un x 


un x 


currying 


cur t 


cur t 


CUT t 


uncurrying 


unc t 


unc t 


unc t 


unique element selection 


theta s 


theta s 


9 s 


element selection 


(added) 


epsilon t 


€ t 


cardinality 


size t 


size t 


size t 


structure 


str t 


(deleted) 


(deleted) 


transitive closure 


t sup -H 


t sup Hr 


/ + 


empty set 


empty 


empty 


0 



TABLE 2. Nonprimitive Extensional Operations: Group 1 



Name 


Old Input Form 


New Input Form 


Publication Form 


pair list 


(x, y) 


(x , y) 


(*, y ) 


left pair section 


(X,) 


(deleted) 


(deleted) 


right pair section 


(>y) 


(deleted) 


(deleted) 


duplication 


DELTA x 


DELTA x 


A X 


membership 


x member t 


x member t 


x e t 


nonmembership 


x nomem t 


x nomem t 


X $ t 


improper subset 


s !subset t 


s .'subset t 


s C / 


proper subset 


s subset t 


s subset t 


s C t 


equality 


s — t 


s — t 


s = t 

cnv G 


converse 


inv t, t sup -1 


cnv t, t sup -1 


domain 


dom t 


dom t 


dom t 


range 


rng t 


rng t 


rng / 


members 


mem t 


mem t 


mem t 


left member 


Lm (x,t) 


x Lm t 


x Lm t 


right member 


Rm (x,t) 


x Rm t 


x Rm t 


member 


Mm (x,t) 


x Mm t 


x Mm t 


right univalent 


run t 


run t 


run t 


left univalent 


lun t 


lun t 


lun t 


bi-univalent 


bun t 


bun t 


bun t 


initial members 


init t 


init t 


init t 


terminal members 


term t 


term t 


term t 


reflexive transitive closure 


t sup * 


t sup ** 


t ' 


domain restriction 


p -> t 


p -> t 


P -> t 


range restriction 


t <- p 


t <- p 


t - p 


restriction 


t restr p 


t restr p 


* ' P 


sequence filtering 


(added) 


p xi t 


p 
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TABLE 3. Nonprimitive Extensional Operations: Group 2 



Name 


Old Input Form 


Nevs Input Form 


Publication Form 


first member 


alpha t 


alpha t 


a t 


last member 


omega t 


omega t 


jj t 


initial sequence 


ALPHA t 


ALPHA t 


A t 


final sequence 


OMEGA t 


OMEGA t 


n t 


ordered union 


t ; u 


t ; u 


t ; u 


cons left 


x cl t 


x cl t 


x cl t 


cons right 


t cr x 


t cr x 


t cr x 


minimum 


min s 


min s 


min s 


maximum 


max s 


max s 


max s 


intersection 


s cap t 


s cap t 


s n t 


set difference 


s \ t 


s \ t 


s \ t 


apply functional record 


t @ hat x 


t @hat x 


t @ X 


apply functional structure 


t ! x 


t ! x 


t ! x 


minimize 


mu t 


mu t 


V t 


database index 


index x d 


x index d 


x index d 


database select 


select x 


x select d 


x select d 


database join 


join x 


x join dblist 


x join dblist 


array to sequence 


as t 


as t 


as t 


sequence to array 


sa t 


t sa i 


t sa i 


seq. to zero-origin array 


saO t 


(deleted) 


(deleted) 


relative product 


rp f t 


t rp f 


t \ f 


relative product inverse 


rpi f t 


f rpi t 


/I t 


array concatenation 


t cat u 


t cat u 


t cat u 


relation sort 


rsort s 


rsort s 


rsort s 


sort 


sort s 


sort s 


sort s 


unit image 


unimg t x 


t unimg x 


t unimg x 


all 


all t 


all t 


all t 


sequence to matrix 


ssm t 


ssm t 


ssm t 



TABLE 4. Primitive Intensional Operations 



Name 


Old Input Form 


New Input Form 


Publication Form 


application 


f @ X 


f @ X 


f @ X 


image 


img f s 


f img s 


i 


f img s 


composition 


f • g 


fog 


t 


f ° 9 


infix to prefix 


(added) 


(op +), (op times), ... 




••• 


left section 


( x +)> ( x -). - 


(lsec x +), (lsec x -), ... 




z+], X- , • • • 


right section 


(+y)» (-y)> - 


(rsec + y), (rsec - y), ... 




+y]» f-y]> • • • 


paralleling 


f 1 g 


f| 1 g 


4 


f II 9 


isomorphism 


f $ t 


f $ t 


f $ t 


formal application 


f @ bar g 


(deleted) 


(deleted) 


functional condition 


(p -> f; g) 


(if p -> f ; g) 


Ip - /; s) 


curry 


curry f 


curry f 


curry / 


uncurry 


uncurry f 


uncurry f 


un curry / 


filtering 


PHI p (d, r) 


p PHI S 


j 


9 $ 5 


iteration 


iter [p -> fj 


(iter p -> f) 


iter p f 


formalization 


4- bar, times bar, ... 


(-*- bar), (times bar), ... 




-F, x , ... 


identity 


Id 


I 


i 



- 19 - 



TABLE 5. Nonprimitive Intensional Operations 



Name 


Old Input Form 


New Input Form 


Publication Form 


while loop 
array reduction 
repeated composition 
value of node 
operate on form 
operate on data 
image of structure 
extension 
restriction 
formal negation 


while [p, f 
f red i 
f sup n 
upsilon f 
phi f 
delta f 
PIf 

extend (t, f) 
restrict (s, f) 
wig p 


(f while p) 
f red x 
f sup n 
upsilon f 
phi f 
delta f 
PIf 

t extend f 
s restrict f 
wig p 


/ while p 

/ § * 
r 

vf 

<pf 

Sf 

n / 

t extend / 
s restrict / 

~ p 



TABLE 6. Miscellaneous Operations 



Name 


Old Input Form 


New Input Form 


Publication Form 


sum 


x + y 


x + y 


X -r y 


difference 


x - y 


x - y 


X - y 


product 


x times y 


x times y 


x x y 


quotient 


x divide y 


x divide y 


x / y 


inequality 


x != y 


x != y 


x # y 


less 


x < y 


x < y 


x < y 


greater 


x > y 


x > y 


x > y 


less or equal 


x <= y 


x <= V 


x ^ y 


greater or equal 


x >= y 


x >= y 


X Z y 


conjunction 


x andsign y 


x andsign y 


xi \ y 


disjunction 


x orsign v 


x orsign y 


x\J y 


negation 


not x 


not x 


—x 


cartesian product 


s cart t 


s cart t 


S X t 



TABLE 7. Data Input Operations and Syntax 



Name 


Input Form 


Publication Form 


identifiers 

strings 

booleans 

relation 

set 

sequence 

list 

subrange set 
subrange sequence 
subrange list 


a, b\ total, etc. 
“abed” 
true, false 
(rel (x : y), ... ) 
(set x y ... ) 

(seq x y ... ) 

(list x y ... ) 
(setrange m to n) 
(seqrange m to n) 
(listrange m to n) 


a. 6 ' , total, etc. 
“abed” 

true, false 

((* *),■■•) 
{x,y, • • • } 

(z, y, ' ) 

<x, y, • • * > 
{m, . . . ,n} 

( m , . . . ,n) 

< m , . . . , n > 
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TABLE 8. RPL Command Types 



Name 


Input Form 


Publication Form 


data definition 


x == y 


* y 


prefix function definition 


f x == y 


} * = y 


infix function definition 


x f y == z 


^ f y = ~ 


write data to a file 


file "name" == x 


file ‘‘name” = x 


read data from a file 


x == (file "name") 


x = file ‘‘name” 


output, form 1 


display x 


display x 


output, form 2 


dis x 


display x 


output, form 3 


d x 


d x 


output, form 4 


X 


X 


output value of definition 


val x 


val x 


output function environment 


env f 


env / 


output entire environment 


env 


env 
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